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INTRODUCTION 


The  diffusion  limited  chemical  equilibrium  combustion  assumption  is  an  important  limit 
case  for  computing  the  chemical  attack  of  reacting  gases  with  a  solid  wall.  Stated  simply,  this 
assumption  means  that  if  the  reacting  gas  can  get  to  the  wall,  it  will  react  completely  with  the 
wall  material  to  produce  equilibrium  products.  If  the  further  assumption  of  complete  similarity 
between  heat  and  mass  transfer  is  made,  i.e.,  unity  Lewis  number,  then  one  of  the  driving 
potentials  for  heat  transfer,  the  hot  wall/gas  phase  enthalpy  can  also  be  computed  using  the 
results  of  a  chemical  equilibrium  assumption.  Taking  the  wall  material  as  fuel  and  the 
combustion  products  of  the  gun  propellant  as  oxidizer,  then  for  a  given  0/F  we  would  like  to 
know  the  temperature  at  which  all  of  the  fuel  is  consumed.  This  temperature-O/F  pair  is  often 
called  the  saturation  point  and  gives  us  the  blowing  parameter  Sq.  The  parameters,  Ba  and  hgw, 
are  used  in  the  following  equations  to  define  the  mass  erosion  rate,  nig. 

PeUgCi^  =q  /(h-h  J 

®®"o  cwrgw  Q'j 

where  Ca„  and  Ch  are  the  Stanton  numbers  for  mass  and  heat  transfer. 

B  o 

Software  and  Engineering  Associates,  Inc.  (SEA)  of  Carson  City,  NV,  has  included  the 
calculation  of  Ba  and  hgw  in  their  version  of  the  NASA  LeRC  CET  (ref  1),  code  which  is  part  of 
their  Material  Ablation  Conduction  Erosion  (MACE)  (ref  2)  set  of  computer  programs.  The 
MACE  code  was  developed  for  reentry  and  rocket  motor  ablation  problems.  In  order  to  extend 
the  code  to  gun  tube  erosion,  the  perfect  gas  assumption  of  the  CET  code  had  to  be  removed  and 
a  better  equation-of-state  implemented  in  the  code.  The  BLAKE  code  (ref  3)  was  subsequently 
used  in  the  gun  community  to  calculate  the  equilibrium  combustion  products  for  gun  propellants 
at  high  pressures,  and  there  are  several  equations-of-state  implemented  in  that  code.  The  choice 
was  then  whether  to  put  the  blowing  calculation  from  our  version  of  CET  into  BLAKE  or  to  put 
the  equations-of-state  from  BLAKE  into  CET.  Because  of  CET's  more  general  chemistry  set,  it 
was  decided  to  add  the  new  equations-of-state  to  CET. 

The  BKW,  Comer/Lennard-Jones  Potential  (LJP),  and  NBS  equations-of-state  were 
added  to  the  MACE/CET  codes  and  checked  out  by  comparing  the  results  with  BLAKE 
calculations  for  the  same  chemical  systems.  Results  from  this  comparison  are  presented  in  the 
verification  section  of  this  report.  The  effects  of  compressibility  on  Ba  and  hgw  for  several 
propellant  systems  and  wall  materials  are  presented  in  the  results  section. 


(2) 

(3) 
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APPROACH 


Equations-of-State 

The  general  equation-of-state  used  here  is  the  same  as  used  in  the  BLAKE  and  Tiger 
codes  (ref  4),  that  is: 


P  =  pZRT  (4) 

where  Z  is  the  compressibility  factor  (1.0  for  an  ideal  gas). 

We  shall  limit  our  brief  discussion  of  the  calculation  of  Z  to  the  Comer/LJP  model,  since 
that  is  the  model  most  used  from  the  BLAKE  code.  References  3  through  5  give  more  complete 
discussions  of  the  various  compressibility  models. 

The  Comer/LJP  model  uses  a  tmncated  virial  equation-of-state,  that  is 

Z  =  [l+pB(T)+/fC(T)]  (5) 

However,  the  Comer  model  assumes  that  the  second  virial  coefficient,  C,  is  a  constant 
and  not  a  function  of  temperature.  The  details  of  the  calculation  of  B  and  C  are  described  in  the 
BLAKE  report  and  are  not  repeated  here.  The  computation  of  the  virial  coefficients  requires  the 
Lennard-Jones  parameters  p  and  CJ;  for  all  the  species  considered  in  the  calculation.  We 
combined  the  Lennard-Jones  parameter  data  from  BLAKE,  CET,  and  SPP  (ref  5)  to  make  as  an 
extensive  a  set  of  data  as  possible.  In  the  case  of  differences  in  the  data,  we  took  the  values  from 
BLAKE  first,  then  CET,  and  lastly  SPP.  For  molecules  where  there  are  no  data,  the  values  for 
nitrogen  are  assumed.  This  last  assumption  is  the  same  as  is  used  in  BLAKE.  Table  1  shows  the 
species  that  are  incorporated  in  the  Compressible  Chemical  Equilibrium  with  Transport 
properties  (CCET)  code. 
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Table  1.  Molecules  with  Lennard-Jones  Parameters  Considered  in  CCET 


AL 

ALCL 

ALCL3 

ALF 

ALF3 

ALN 

ALO 

ALS 

AL2 

AIR 

AR 

ASH3 

B 

BBR3 

BCL 

BCL2 

BCL3 

BF 

BF2 

BF3 

BIS 

BO 

B(OCH3)3 

B2 

B2H6 

B203 

BE 

BEBR2 

BECL 

BECL2 

BEF 

BEF2 

BEE 

BE2 

BR 

BRF 

BRF3 

BRO 

BR2 

C 

CBRF3 

CBR4 

CCL 

CCLF3 

CCL2 

wimssmmm 

CCL3 

CCL3F 

CCL4 

CF 

CF3 

CF4 

CH 

CHBRCLF 

CHBRCL2 

CHBR3 

CHCLF2 

CHCL3 

CHF3 

CH2BRCL 

CH2CLF 

CH2CL2 

CH2F2 

CH2E 

CH3BR 

CH3CL 

CH3F 

CH3I 

CH30H 

CH4 

CN 

CO 

COS 

C02 

CP 

cs 

CS2 

C2 

C2H2 

C2H4 

C2H6 

C2H5CL 

C2H50H 

C2N2 

CH30CH3 

CH2CHCH3 

CH3CCH 

CYCLO-C3H6 

C3H8 

N-C3H70H 

CH3COCH3 

CH3COOCH3 

N-C4H10 

ISO-C4H10 

C2H50C2H5 

CH3COOC2H5 

N-C5H12 

C(CH3)4 

C6H6 

C6H12 

N-C6H14 

CD 

CL 

CLCN 

CLF 

CLF3 

CLO 

CL2 

F 

FCN 

F2 

H 

HBR 

HCN 

HCL 

HF 

HI 

HS 

H2 

H20 

H202 

H2S 

HE 

HG 

HGBR2 

HGCL2 

HGI2 

I 

ICL 

12 

KR 

LI 

LBR 

LICN 

LICL 

LIF 

LE 

LIO 

LI2 

LI20 

MG 

MGCL 

MGCL2 

MGF 

MGF2 

MG2 

N 

NFS 

NH 

NH3 

NO 

NOCL 

N2 

N20 

NA 

NABR 

NACN 

NACL 

NAF 

NAI 

NAO 

NAOH 

NA2 

NA20 

NE 

0 

OF 

0F2 

OH 

02 

P 

PCL 

PCL3 

PF 

PF3 

PH3 

PN 

PO 

PS 

P2 

P4 

S 

SF6 

SO 

S02 

S2 

S2F2 

SI 

SICL 

SICL4 

SIF 

SffCLS 

SIF2CL2 

SIF3CL 

SIF4 

SIH4 

SIO 

SI02 

SIS 

SI2 

SNBR4 

SNCL4 

UF6 

XE 

ZN 
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Verification 


The  method  of  verification  of  the  modifications  made  to  the  CCET  code  was  to  compare 
the  new  code  with  results  from  both  the  BLAKE  and  CET  codes.  Since  the  compressibility 
models  were  taken  directly  from  the  BLAKE  code,  the  results  should  compare  directly  with  that 
code.  Comparisons  with  the  CET  code  would  establish  that  we  had  not  broken  the  CCET  code 
when  we  implemented  the  compressibility  computations  and  would  establish  that  all  three  codes 
computed  the  same  result  when  compressibility  was  not  a  factor. 

Figure  1  shows  a  comparison  of  selected  species  mole  fractions  as  predicted  by  BLAKE 
and  CCET  for  a  pressure  of  5000  atmospheres  and  a  temperature  range  from  1000  to  5000°K. 
The  upper  set  of  symbols  refers  to  the  BLAKE  numbers,  and  the  lower  set  to  the  CCET  numbers. 
Figures  2  and  3  show  a  comparison  for  compressibility  and  enthalpy.  In  Figure  2  the  double 
valued  point  for  the  BLAKE  compressibility  stems  from  an  apparent  convergence  tolerance 
problem  when  the  results  are  computed  going  up  a  temperature  isoline  as  compared  to  going 
down  a  temperature  isoline.  The  input  file  listings  for  the  BLAKE  and  CCET  data  are  shown  in 
Tables  2  and  3,  respectively. 

M30A1 


CO 

N2 

H20 

. -X . . 

H2 

C02 

CH4 

— ^ — 

OH 

CH3 

H 

0 

CO 

N2 

. 

H20 

••  . 

H2 

. . 

C02 

CH4 

OH 

CH3 

H 

0 

Figure  1 .  Comparison  of  BLAKE  and  CCET  species  mole  fractions  for  M30A1  propellant. 
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M30A1  Compressibility 


— CCET 
-■—BLAKE  UP 
BLAKE  DOWN 


T(K) 


Figure  2.  Comparison  of  BLAKE  and  CCET  compressibility  for  M30A1  propellant. 


M  30 A  1 


T  (K) 

Figure  3.  Comparison  of  BLAKE  and  CCET  gas  phase  enthalpy  for  M30A1  propellant. 
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Table  2.  Input  File  for  CCET/M30A1  Propellant 


#  Test  of  CCET  versus  Blake 

#  a)  Assigned  temperature  and  pressure  problem  (tp  =  t). 

#  b)  Pressure  is  5000  atm  (atm  =  t). 

#  c)  Temperatures  are  given  in  degree  K  (t  =  ...)• 

#  d)  M30A1  from  LIBRARY  cards. 

# 

$DATA 

LIB  =  T,  MOLE  =  F,  LEBNAM  =  'NCI 260',  'NG',  'NQ',  'EC,  'KS' 
MASMOL  =  27.9,  22.42,  46.84,  1.49,  1.0,  .25,  .1, 

SEND 

'ALC,  'C, 

OMIT 

02 

C(GR) 

C2N 

C2H 

OMIT 

C2 

CH20 

N02 

H2S 

OMIT 

S20 

S02 

K(S) 

KOH(S) 

OMIT 

K20 

K202 

K02 

OMIT 

C 

CH 

K2 

OMIT 

KCO(S) 

KSO(S) 

K20(S) 

NA2(S) 

OMIT 

C2H4 

C2N2 

C2H2 

CH2 

OMIT 

HNO 

HN03 

K2S(S) 

OMIT 

NAMELISTS 

$CET 

LJP  =  T, 

TP  =  T, 

T  =  1000,  2000,  3000,  4000,  5000 

P  =  5000,  ATM  =  T, 

SEND 
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Table  3.  Input  File  for  BLAKE/M30A1  Propellant 


TIT,  M30A1  PROPELLANT  UP 
REJ,  02,  C(S) 

REJ,  C2N,  C2H,  C2,  CH20,  N02,  H2S,  S20,  S02,  K$,  KOH$,  K20, 
REJ,  K202,  K02,  HN03,  C,  CH,  K2,  N, 

REJ,  KCO$,  KSO$,  K20$,  NA2$ 

REJ,  C2H4,  C2N2,  C2H2,  CH2 
REJ,  HNO,  HN03 
REJ,K2S$ 

0RD,N2,C0,H20,K0H,HS 

RET 

CM2,  NCI 260,  27.90,  NG,  22.42,  NQ,  46.84,  EC,  1.49, 

KS,  1.0,  ALC,  .25,  C,.l 
ISOline,  P,  5000,  T,  1000., 4,  5000. 

STOP 


The  above  comparison  was  repeated  for  M919  propellant  for  a  pressure  of  3000 
atmospheres  and  temperatures  from  1000  to  4000°K.  The  results  were  comparable  to  the 
M30A1  comparisons. 

The  final  comparison  was  at  one  atmosphere  using  a  stoichiometric  mixture  of  ethanol 
and  oxygen  at  a  temperature  of  3000°K.  Comparisons  between  BLAKE,  CCET,  and  CET  were 
essentially  identical. 

RESULTS  AND  DISCUSSION 

The  purpose  of  this  work  was  to  develop  an  equilibrium  code  that  would  compute  values 
of  Ba  and  hgw  for  use  in  thermochemical  erosion  studies,  which  include  the  effects  of  nonideal 
gas  equations-of-state,  that  is,  compressibility.  Results  were  computed  for  the  M30A1  propellant 
in  conjunction  with  chromium  with  the  compressibility  option  of  CCET  turned  on  and  off  for  a 
pressure  of  3000  atmospheres.  The  values  for  Ba  were  found  to  be  essentially  identical  for  the 
two  cases  and  hence  independent  of  compressibility.  These  results  for  the  gas/wall  enthalpy,  hgw, 
were  significantly  different  and  are  shown  in  Figure  4. 
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These  results  show  that  the  effects  of  compressibility  are  important  when  computing  the 
heating  and  diffusion  rates  for  thermochemical  ablation/erosion  for  gun  tube  applications. 

SUMMARY  AND  CONCLUSIONS 

Modifications  to  extend  the  NASA  LeRC  GET  to  include  compressibility  effects  have 
been  demonstrated.  The  new  code  has  been  named  the  Compressible  Chemical  Equilibrium  with 
Transport  properties  (CCET)  code.  Good  agreement  between  the  gun  community's  BLAKE  code 
and  the  CCET  has  been  shown. 

The  effect  of  compressibility  on  thermochemical  erosion  potentials,  Ba  and  hgw,  has  been 
calculated  for  two  propellants,  M30A1  and  M919.  It  was  concluded  that  compressibility  is 
important  in  computing  thermochemical  erosion  rates  in  gun  tubes. 

The  compressibility  models  in  BLAKE  have  been  found  to  be  adequate  for  internal 
ballistics  computations.  However,  for  other  uses,  we  feel  that  more  accurate  models  should  be 
incorporated  into  the  CCET.  The  calculation  of  diffusion  coefficients  and  Schmit  numbers  for 
computed  equilibrium  conditions  is  also  important  for  many  applications  and  should  be  added  to 
the  CCET  code. 
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AMSTA-AR-AET-0,  BLDG.  183  1 

AMSTA-AR-FSA,  BLDG.  61  1 

AMSTA-AR-FSX  1 

AMSTA-AR-FSA-M,  BLDG.  61  SO  1 

AMSTA-AR-WEL-TL,  BLDG.  59  2 

PICATINNY  ARSENAL,  NJ  07806-5000 

DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-DD-T,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MD 
21005-5066 


DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-WM-MB  (DR.  B .  BURNS)  1 

ABERDEEN  PROVING  GROUND,  MD 
21005-5066 

COMMANDER 

U.S.  ARMY  RESEARCH  OFFICE 

ATTN:  TECHNICAL  LIBRARIAN  1 

P.O.  BOX  12211 

4300  S.  MIAMI  BOULEVARD 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SIORI-SEM-L  1 

ROCK  ISLAND,  IL  61299-5001 

COMMANDER 

U.S.  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIBRARY)  1 

WARREN,  MI  48397-5000 

COMMANDER 

U.S.  MILITARY  ACADEMY 

ATTN:  DEPT  OF  CIVIL  &  MECH  ENGR  1 

WEST  POINT,  NY  10966-1792 

U.S.  ARMY  AVIATION  AND  MISSILE  COM 
REDSTONE  SCIENTIFIC  INFO  CENTER  2 

ATTN:  AMSAM-RD-OB-R  (DOCUMENTS) 
REDSTONE  ARSENAL,  AL  35898-5000 

COMMANDER 

U.S.  ARMY  FOREIGN  SCI  &  TECH  CENTER 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT.  AND  ENGINEERING  CENTER, 
BENET  LABORATORIES,  CCAC,  U.S.  ARMY  TANK-AUTOMOTIVE  AND  ARMAMENTS  COMMAND, 

AMSTA-AR-CCB-0.  WATERVLIET,  NY  12189-4050  OF  ADDRESS  CHANGES. _ 


